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ABSTRACT 


Image differencing has been investigated as a technique for use with Landsat 
digital data to delineate areas of land cover change in an arban environment. 
Landsat data collected in April 1973 and April 1975 for Austin. Texas, were 
geometrically corrected and precisely registered to United States Geological 
Survey (USGS) 7.5-minute quadrangle maps. At each pixel location reflectance 
values for the correspmiding bands were subtracted to produce four difference 
images. Areas of major reflectance differences are isolated by thresholding 
each of the difference images. The resulting images are combined to obtain 
an image data set of total change. These areas of reflectance differences 
were fotind, in gemral, to correspond to areas of land cover change. Informa- 
tion on areas of land cover change was incorporated into a procedure to m?sk 
out all nonche?ge sreas and perform an unsupervised classification only for 
data in the change areas. This procedure identified three broad categories: 

(1) areas of high reflectance (construction or extractive). (2) changes in agri- 
cultural areas, and (3) areas of confusicm between agricultural and other areas. 
Efforts are underway to include 1977 Landsat data and to rigorously evaluate 
the results of the image differencing technique and the unsupervised classifica- 
tion. 
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SECTION 1 - BACKGROUND 


In 1975 a Joint reaearch project between Goddard ^pace Fll^t Center (GSFC) 
and the Bureau of the Cenaus began. The program la to examine the utility of 
Landaat remotely aenaed data for monitoring land cover/land uae activity 
around major urbanized areaa in the United Statea. Early work in the project 
has shown that major land cover classes in the urban fringe (transition zone be- 
tween urban and rural populations) can be identified from satellite data (Refer- 
ence 1). 

Land cover change is a key element in the uriMtnized area boundary delineation 
program conducted by the Census Bureau. The frequent* synoptic nature of 
Landsat data makes it well suited for the detection of land cover change. In an 
extension of the earlier work on this project* a technique for the delineation of 
areas that have undergone major changes in spectral reflectance (indicating 
changes in land cover) has been developed. 

This document reports progress to date in the development of one af^roach for 
monitoring land cover change and also describes objectives of future research. 
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SECTION 2 - ANALYSIS TECHNIQUE 

The fundamental operatira in automated change detection is the subtraction of 
one digital image from another. The differencing (^ration is carried out by 
subtracting the reflectance value for a point in one image from the reflectance 
value for the same location in the second image. Landsat I or II multispectral 
scanner (MSS) output consists of four bands or digital Images. Thus, to differ- 
ence all of the Landsat MSS data in a scene requires four operations per pixel 
location— one for each of the four MSS bands. 

The subtracti<« of one digital image from the other results in an image in which 
positive and negative values represent areas of change and zero values indicate 
areas that remained unchanged. When applied to Landsat data with a data range 
of 0 to 127, the potential range of difference values is -127 to -•-127. In the 
Small htteractive Image Processing System/Video Image Communication and 
Retrieval System (SMIPS/VICAR) system (Reference 2), all results are stored 
as positive values scaled from 0 to 255 by the addition of a constant. This oper- 
ation can be expressed mathematically as 

^ijk ^ 


where AX = change image 
X(l) s image at time 1 
X(2) = image at time 2 
C > a constant 
i a 1. . . n number of lines 
j = 1. . . n number of columns 
k s 1. . . n number of bands 

When applied to Landsat data, this results in four different images that can be 
subsequently analyzed using standard image processing techniques. 


For this project areas of change were identified for each of the four MSS bands 
and these results combirod to obtain an image representing total change. 

Change areas for each band of the difference im^e were determined by identify- 
ing those pixels with intensity values above or below calculated threshold levels 
(i.e. , density slicing the image). Pixels in the difference image with reflec- 
tances above the high threshold represent areas that have increased in reflec- 
tance from time 1 to time 2 (high reflectance change), whereas pixels with 
values below the lower threshold represent areas that have decreased in reflec- 
tance from time 1 to time 2 (low reflectance change). Pixels with reflectances 
that fall between the threshold values are considered to have remained un- 
changed. This procedure results in a series of images depicting areas of 
change for each band which can then be combined to obtain an image showing all 
areas of change. 

This process makes two assumptions: (1) changes in land cover Aand use will 
result in corresponding changes in reflectance values, and (2) these reflectance 
changes will be of sufficient magnitude to distinguish them from the normal var- 
iation of features within the images. 
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SECTION 3 - STUDY SITE 


The test site for tiiis stucty is an area of approximately 1200 square kilometers 
encompassing Austin, Texas, and its envU-ons. This area was chosen because 
the Census Bureau conducted pretests there for the 1980 census at the start of 
this woric. The Bureau intended to compare pretest results with satellite data 
products. As part of that earlier work, 1973 and 1975 images were registered 
to United States Geological Survey (USGS) topograi^ic maps. The availability 
of registered images and the interest of the Bureau of tte Census led to the use 
of Austin as the site for image differencing/change detection eiqierimentation. 
To date, however, results of the census pretest have not been made available. 

Austin is located in central Texas between the Coastal Plain and the Edwards 
Plateau. Its terrain varies from level, c^n farmland and fields in the east to 
tree*covered hills in the west. The Austin central busii^ss district (CBO) is 
suzTounded by commercial, industrial, and high>density residential develop- 
ment. Further from the center city—primarify to the north, east, ami south- 
lie low-density residential development and agricultural areas. In addition, a 
limited number of estates and trailer parks are mixed into this diverse land- 
scape. 
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SECTION 4 - DATA 


Subsets of Landsat data collected on March 17, 1973 (Figxire 4-l(A)), and 
March 24, 1975 (Figure 4-l(B)), were used as the principal input for this study. 
Near-anniversary data sets were selected to minimize phenological and sea- 
sonal variations. Each subimage was geometrically corrected and precisely 
registered to USGS 7-1/2-minute quadran^e maps and resampled to 60-metei- 
square picture elements (Figures 4-2(A) and (B)) using the Digital Image Recti- 
fication System (DIRS) implemented at GSFC (Reference 3). The SMIPS/VICAR 
system available at GSFC was used for all subsequent processing and analysis. 

Supporting data for the study included USGS topographic maps (photo-revised in 
1973), low-altitude black and white photography taken in 1973 and 1974, a few 
frames of low -altitude black and white photography taken in 1975, and 1977 high- 
altitude black and white imagery. 
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MARCH 17, 1973 


MARCH 24.1975 


Figure 4-1. 


La.asat scenes U37-1637S and 7374-16,30. me AusUn. Texas. 
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Figure 4 -2. Enhanced Pictures of the Geometrically Corrected and Registered 
Portions of Landsat Scenes 1237-16372 and 1974-16130 
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SECTION S - RESULTS AND DISCUSSION 

For this itucly a dlfforonoo image waa produced by eid>tractiag 1973 pixel valuea 
from the oorre^oodiog elemeata ia the 1975 image and adding a o<mstant of 
128. Thua» Equation (2-1) beoomea 

The four remlting difference images are the baele of all further analysis. 

5. 1 DIFFERENCE DATA STATISTICS 

Differencing procedure results are illustrated by histograms of MSS 7 for the 
1973* 1975, and difference data eets (Figure 5-1). The histogram of the differ- 
enced data (Figure 5-l(C)) is centered near 128 with a much tighter standaid 
deviation than the histograms of the raw data. The histograms of all four dif- 
ference images (Figure 5-2) are ti^tly clustered near 128. The means and 
standard deviations for each of these images and their raw data counterparts 
are presented In Table 5-1. Mean values near 128 are expected for the differ- 
ence data sets because most areas have remained unchanged during the 2-year 
period. 


Table 5-1. Means and Standard Deviations of 1973, 1975, 
and Differenced Images 
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Figure 5-1. Histograms of MSS 7 for 1973, 1975, and the Difference Data Set, 
With Mean Values Shown for Each Distribution 
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Figure 5-2. Histograms of the Four Difference Data Sets, With Mean Values Shov n 


Mean values and standard deviations for MSS-6 and MSS-7 difference data are 
hig^r and broader, reqsectively, than the means or standard deviations in 
MSS-4 and MSS-5 difference Images. This suggests that an overall shift in ra- 
diance has occurred for these spectral bands. Indeed, rainfall datn for March 
1973 and March 1975 (Figure 5-3) show a cumulative rainfall of 1.2 inches for 
the 2 days prior to the March 17, 1973, image, whereas no ai^reciable rainfall 
was recorded for the 10 days prior to the March 24, 1975, image. Such a dif- 
ference In ground moisture would contribute to the reflectance difference be- 
tween the images, particularly for MSS 6 and MSS 7, because both bare soil 
and healthy vegetati(» decrease their reflectance of infrared radiation as their 
moisture content increases. Impervious areas, like the CBD, would not exhibit 
such changes. Given the difference in rainfall, resp<mses for MSS 6 and MSS 7 
for 1973 would be expected to be more varied and on the average lower than the 
values for 1975. Indeed, the mean values and standard deviations for MSS (i 
and MSS 7 of the difference image (1975-1973) are higher and broader than 
MSS 4 and MSS 5 of the same product. 

5. 2 DIFFERENCE THRESHOLD 

A direct display of the difference data produces a flat, low-contrast picture that 
is difficult to interpret (Figure 5-4). The intensity of an area in the image indi- 
cates the magnitude of the reflectance difference between the two dates for that 
location. Areas that have ejq)erienced major change appear light (increased re- 
flectance) and dark (decreased reflectance) in the image. These light and dark 
areas correspond to the data located in the tails of the difference data histo- 
grams. The pixel locations corresponding to these differences can be selec- 
tively displayed by density slicing each of the difference images. 

In this experiment the threshold values for the density slice were set at three 
standard deviations above and below the mean difference value (Table 5-2). 
Initially, threshold values of plus and minus two standard deviations were se- 
lected, but this appeared to overestimate the probable amount of change. 
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Figure 5-3. Daily Rainfall Data for Austin, Texas, 
During March 1973 and March 1975 




Figure 5-4. 1975-1973 Difference Images by Band 
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Table 5-2. Threshold Values for Density SliM 


DiFFERENC£D 
MSS BAND 

3 STANDARD 
DEVIATIC^ 

-3 STANDARD 
DEVIATIONS 
(LOST THRESHOLD) 

*3 STANOASO 
OEVIATIOUtt 
(HIGH THRESHOLD) 

4 

ao.i 

106 

149 

5 

27.9 

101 

157 

S 

30.3 

101 

162 

7 

36.7 

95 

166 


Picture elements with differences less than, or greater than, the threshold 
value are displayed in Figure 5-5. The top row of images (Figures 5(A) through 
(D)) shows those pixels that have decreased in reflectance from 1973 to 1975, 
for each of the MSS haras. The bottom row of images (Figures 5-5(E) through 
(H)) displays those pixels that have increased in reflectance. The areas delin- 
eated in these images correspond to all locations diat had a variation in meas- 
ured reflectance greater than the noise threshold of three standard deviations. 
These images were combined to generate a series of data sets that depict areas 
of decreased reflectance (low change), increased reflectance (high change), and 
decreased and increased reflectance (total change) (Figure 5-6). 

5.3 PRELIMINARY EVALUATION 

The areas delineated in this way may provide useful information for identifying 
certain types of land cover change. Preliminary comparisons indicate that 
areas which have decreased in reflectance (low reflectance change) are primar- 
ily vegetated or agricultural areas. A^eas of increased reflectance, however, 
seem to include both agricultural areas and areas of construction activity. A 
comparison of the low-reflectance change image and high-reflectance change 
image shows that 407 pixels increased in reflectance in one or more bands and 
also decreased in reflectance in one or more of the remaining bands. A color 
display of the total change image (Figure 5-7) highlights these areas as white, 
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Figure 5-5. Areas of Reflectance Difference Greater Than Three Standard Deviations 
Above and Below Mean Values by Band 
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Figure 5-6. Composite Images ol cnange Areas 
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Figure 5-7. Total Change for the Scene. Low change is shown in red, 
high change in blue. Areas with both low and high change 
in various bands are shown in white. 
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with low oluu^ <Ufptoyed in red and high change dtqdayed tn blue. Areas of 
both low and high reflectance change were, In teot, limited to agricultural areas 
and mi^ be die result of a unique combination of agricultural practices and en- 
vironmental conditions. 

The image of total change was combined with a color enhancement of the raw 
data to produce an image that facilitates the evaluation of the change areas by 
providing a hrame of reference for these areas (Figure 5-8). Preliminary com- 
parisons of the change areas with supporting information are encouraging. Sev- 
eral major features have been correctly identified as change. Many of the 
change pixels are associated with agricultural areas, as tyjdfled by the agricul- 
tural areas adjacent to the river in the eastern portitm of the image. Other 
major chai^ areas are the result of construction activity, examples of which 
are freeway constructlmi and prdOable residential building (butterfly shape) in 
the northwestern portion of the study area. 

At present only general comparisons can be made regarding the accuracy of the 
change areas due to the lack of coincidental ground-truth information. These 
preliminary comparisons serve to demonstrate the application of such a proce- 
dure as a first-order estimate of land cover change. The principal advantage 
of this approach is that, given registered Imagery, it is a computationally sim- 
ple technique that can be accomplished with a minimum of analyst effort. 

5.4 CHANGE AREA CLASSIFICATION 

Thus far the procedure has identified areas of change within the image but has 
not provided specific information regarding the typ>es of change that have 
occurred. Given the first-order estimate of land cover change, this informa- 
tion has been Integrated Into a procedure for classifying tttose diange areas. 

The change data set can be used as a binary mask to extract the data for areas 
of change and thereby reduce the amount of data that must be processed by 
97 percent. 
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Total Change Displayed Against a Background of the 1975 Data 
Total change is shown In yellow. 
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^oh a change area ma^ e^a applied to the 1976 <btta and to the differmce data 
seta. The resultli^ data sets, omidsting of response values for the areas of 
change and aero values elsewhere, were then ana^zed usii« the SMIPS/VICAR 
version of die HINDU clustering algorithm. The q>tlon exists within the pro- 
gram to exclude specified values from the calculations. Utilizing this function 
to eliminate zero-valued pixels, processing was only carried out on the data 
for chan^ areas. This an>roach has the advantage of substantially reducing 
the processing time. Further, the level of discrimination achieved by cluster- 
ing can be expected to be more detailed for these areas tl.an if all the data for 
the test site were processed. The initial classiflcations on both the 1975 data 
and the change area data are encouraging. 

The clustering of the 1975 data in areas of change resulted In seven categories. 
Due to the lack of available ground truth, these results have not undergone de- 
tailed analysis. However, a preliminary evaluation grouped six of the cate- 
gories into two general classes: high-reflectance change areas (primarily 
construction and extractive) (3213 points or 1.3 percent) and agricultural areas 
(3664 points or 1.4 percent). The seventh category combined both agricultural 
and noni^icultural areas (1473 points of 0. 5 percent). These three general 
classes are shown in Figure 5-9. Results of the clustering of the difference 
data have not yet been evaluated. 
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Figure 5-9. Land Cover Classification of 1975 Data for Change Areas. 

Agriculture Is shown In green, highly reflective areas in 
red, and confusion areas In yellow. 
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SECTION 6 - FUTURE WORK 


The limited ground truth for 1975 precluded e ligoroue evaluation of the awu- 
racy of the procedure. Aerial photo coverage is available for 1977. as is a 
clear Landsat image for April 10, 1977. To test more thoroughly the previ- 
ously described technique, it was decided to register the 1977 image to the 
earlier scenes and to rvpeat the change detection process using the 1973 and 
1974 data sets. The 4-year time span of these images approximates the 
S-year interval under which the Bureau of the Census will operate. 

In addition to evaluating change over a 4-year period, the use of the 1977 data 
set provides an qnxx'bmity for a comparison of change over each 2-year period 
as a possible indicator of land cover chat^e patterns. 

6. 1 SCENE REGISTRATION 

Registration of the 1977 image to the earlier data set will be performed using 
a combination of DIRS and VICAR routines. This processing is nearly com- 
pleted. Only the resampling of the 1977 image to conform with the 1973 data re- 
mains to be done. 

DIRS was used to perform a flrst-oider correction for major sources of dis- 
tortion (Earth rotation skew, scanner over-sampling, mirror velocity, and so 
forth) and to resamide the data for 60-meter-square pixels, the resolution of 
the 1973 and 1975 registered images. The first-order correction results in 
data sets that approximate the projection of the registered imagery. 

The near-alignment of the data sets facilitates the identification of correlation 
points to be used for precise alignment of the images. Initially, 70 such points 
were selected. Two points were subsequently eliminated due to poor scene-to- 
scene correlation. The correlation values obtained from registration runs 
yielded a typical residual of 0. 2 pixel between the tie points and the resampling 
grid surface. This value, combined with the use of 68 tie points for an area 
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only 640 lines by 515 samides, suggests that the 1977 image will register verj' 
accurately to the 1973 data. 

6.2 GROUND TRUTH 

Ul^SS 7-1/2-minute topographic m^ps, photo-revised in 1973, serve as the 
basis for comparison with the 1973 image. These maps are being updated 
using 1977 black and white aerial photography to provide a source of informa- 
tion concurrent with the 1977 Landsat data. Together these maps will provide 
additional land cover information that will provide a basis for evaluating the 
satellite-derived land cover change. 

To facilitate the comparison of the photo -interpreted and the satellite-derived 
change areas, the photo-inteipreted change information will be digitized using 
the television scanner input capabilities of the Atmospheric and Oceanc^raphic 
Image Processing System (AOIPS) Hazeltine terminals. The resulting digital 
data will be corrected for scanner distortion, scaled, and registered with the 
satellite data. 

6. 3 DATA COMPARISON 

The availability of registered satellite data and ground-truth information will 
permit an automated point-by-point comparison of the data results. This will 
make it possible to evaluate the results at each step in the differencing se- 
quence (e. g. , evaluation of change by band and evaluation of total change). The 
digitized ground truth will also make it possible to quickly evaluate the utility 
of the : NDU clustering program to classify change areas. Such rigorous analy- 
sis of ti various products is only practical when an automated technique such 
as this is employed. 
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SECTION 7 - SUMMARY 


An efficient image differencing scheme for the detection of major changes in 
land cover has been demonstrated. TIm procedure subtracts the data value at 
a point from the data value for the same location in the second scene, then adds 
a ccHistant to the differem:e. These results can be displayed in photographs 
using existing image processing systems, but their interpretation is not ob- 
vious. To facilitate identification of change areas, file difference images are 
density sliced at three standard deviations above and below their mean values. 
Pixels that show a variation greater than the threshold are taken to be areas 
of land cover change. Hie combination of threshold results for each band 
produces a total change image for the scene. 

Besides providing a map of areas with substantial reflectance changes, the total 
change image can be used as a binary mask to extract raw reflectance data for 
change areas. This procedure reduces the amount of raw data for further analy 
sis by approximately 97 percent. The reduction permits the HINDU clustering 
algorithm to rapidly process the scene. Such a classification of 1975 data re- 
vealed that most areas of high reflectance change are the result of construction 
activity. In the Census application these areas could be considered areas of 
residential development. 

To test the image differencing technique further, a 1977 Landsat scene will be 
registered to the earlier data sets and used to analyze change in the study site 
over a 4-year period. This work has two primary benefits: (1) the 4-year 
time span more closely approximates the 5-year timeframe within which 
the Bureau of the Census will update its records, and (2) good ground truth is 
available for evaluating the results. 

Using 7-1/2-minute quadrangle sheets (photo-revised in 1973) and 1977 aerial 
photos, it will be possible to derive ground -truth land-cover change informa- 
tion. This data will be digitized using a television scanner available on the 
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AOIPS Hazeltine terminals and then registered to the satellite data. Such a dig- 
itized ground-truth data base will permit extensive point-by-point evaluation of 
the Landsat-derived change areas and of the HINDU land cover classification. 

The preliminary results are encouraging. They show that differencing Landsat 
images is an efficient technique to delineate areas of land cover change. Fur- 
ther, the information obtained from such a procedure can be easily integrated 
into an analysis and classiflcaticm procedure to identify land cover chaise by 
type. 
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